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Accreting intermediate-mass black holes (IMBHs) have been proposed as an explanation for ultraluminous X-ray 
sources (ULXs). Runaway core collapse inside a massive cluster is a possible mechanism for IMBH formation. 
But if so, why are ULXs only rarely found associated with a cluster? We use a simple analytical approximation 
to show that rapid core collapse can occur in two physical regimes. For cluster masses ~ 10 6 Mq, an IMBH 
may be formed if the collapse timescale is ^ 3 Myr, as already well known; the parent cluster is expected to 
survive. For cluster masses ~ 10an IMBH may result from a core collapse on even shorter timescales 
(« 0.5 Myr), when the protocluster is still embedded in optically thick gas. Most clusters in this latter regime 
are disrupted “explosively” as soon as the gas is ionized by the OB stars. We speculate that this process may 
produce isolated ULXs with masses up to a few 10 2 Mq, surrounded by a loose OB association, and perhaps by 
a nebula, remnant of the expanding gas from the disrupted protocluster. 


1. INTRODUCTION 

Various alternative models have been proposed to 
explain the nature of ULXs in nearby galaxies. In¬ 
homogeneous accretion flows may allow luminosities 
a factor of 10 in excess of the Eddington limit ( 0 )- 
Moderate geometrical beaming (Q) in the direction of 
the observer may boost the observed brightness by an 
order of magnitude. Thus, sources with an apparent 
isotropic X-ray luminosity up to a few 10 39 erg s~ 4 can 
be explained by stellar-mass black holes (BHs) with 
masses « 5-15 M 0 , without the need to invoke new 
astrophysical mechanisms. However, these scenarios 
become problematic for ULXs with steady apparent 
X-ray luminosities > 10 4 ° erg s~ 4 : they require either 
strong beaming, such as a relativistic jet pointing to¬ 
wards the observer ( 0 ; 0 ) , or a more massive ac- 
cretor, with M ~ 100-1000 M 0 (IMBHs: Q). It is 
possible that ULXs are a mixture of different classes of 
objects. Observational and statistical arguments have 
been made against a predominance of beamed sources 
( 0 ). In other cases, circumstantial evidence in favor 
of IMBHs has been proposed from X-ray spectral and 
timing studies ( 0 )- 

One of the difficulties of the IMBH scenario is how 
to introduce a plausible formation mechanism. It 
has been suggested that IMBHs could be the rem¬ 
nants of very massive 1000 M 0 ) zero-metallicity 
Population-Ill stars ( 0 ). However, this is difficult 
to reconcile with the relative abundance of ULXs in 
young star-forming regions, especially in starburst and 
merging galaxies. Another possibility ([9f) is that at 
least some ULXs may be the nuclei of smaller satel¬ 
lite galaxies accreted and tidally disrupted by a larger 
host (Hq)). 

Alternatively, it has been proposed (|llj, 10, [H|) 
that IMBHs could be formed in the dense cores of 
massive, young stellar clusters or super star clusters. 
Analytical approximations and numerical simulations 
have shown that, for a suitable range of cluster masses 


and densities, mass segregation and the Spitzer insta¬ 
bility (Q) can lead to a runaway core collapse on a 
timescale shorter than the lifetime of its most massive 
stars (ps 3 Myr). The same simulations also show that 
the mass of the collapsing core is ~ 10 -3 times the 
total mass of the cluster. Thus, this scenario could 
explain the formation of IMBHs with a mass up to 
~ 1000 M 0 . 

The main difficulty of theoretical models invoking 
IMBH formation inside a cluster is that, in fact, most 
ULXs do not have a bright, massive cluster as their 
optical counterpart, apart for few notable exceptions 
(e.g., a bright ULX in M 82 is associated with the star 
cluster MGG-11: 0)- I n some cases, such as the 
ULXs in the Cartwheel Galaxy f [lfijf , the distance is 
such that we cannot draw firm conclusions on the sig¬ 
nificance of possible ULX-cluster coincidences. In the 
starbursting Antennae Galaxies, most ULXs are dis¬ 
placed by ss 300 pc from their nearest star clusters 
(Eii M)- I n the starburst galaxy NGC7714, no op¬ 
tical counterparts have been found for the two bright¬ 
est ULXs ([Ts . [u|) : there are no objects brighter 
than My ~ Mb ~ —9 within a few hundred pc of 
the two sources. This rules out young clusters with 
masses 10 4 M 0 . In other cases (e.g., in NGC5408: 
H3; in NGC5204: 0|; in NGC4559: [z|; in Holm- 
berg IX: [23|) the optical counterpart is thought to be 
an individual O or B star. 

One explanation for the lack of ULX-cluster associ¬ 
ations is that the accreting systems are runaway bina¬ 
ries ejected from a cluster. However, the kick veloc¬ 
ities required to explain displacements of a few hun¬ 
dred pc would rule out BH masses 20 M 0 (HU). 
Hence, the ejection scenario is not consistent with 
IMBHs. An alternative possibility is that the par¬ 
ent cluster has already dispersed. Physical mecha¬ 
nisms leading to the expansion and disintegration of 
star clusters are tidal disruption, or mass loss via SN 
explosions and winds from evolved massive stars (e.g., 
H3; @|)- However, it is doubtful that these processes 
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can entirely dissipate a massive cluster on timescales 
as short as ~ 10 ' yr. 

In this paper, we use simple analytical approxima¬ 
tions to compare the timescale for IMBH formation 
with characteristic timescales of cluster evolution. We 
show that for some density profiles, the collapse can 
already occur during the initial protocluster phase, 
when the stars are still embedded in optically-thick 
gas. We then show that for a range of masses and 
densities, the parent clusters may not survive the em¬ 
bedded phase, leading to the formation of an appar¬ 
ently isolated IMBH. 

2. VELOCITIES AND DENSITIES IN A 
CLUSTER CORE 


for the more centrally concentrated Wo = 9 profile. 

We could approximate the Coulomb logarithm 
In A « 7-9 for the parameter range of interest here. 
More accurately, we use the definition of A (|b||, 
chapter 2), to obtain A = 37\af/(Gm) « OAM/m, 
where r h is the half-mass radius. Expressing M oc 
jO - 1 / 2 (0)o-3(0), where the proportionality constant de¬ 
pends on the cluster model, we obtain: 

A = aG _ 3 / 2 p _ 1 / 2 (0)o-?(0)m _ 1 (0). (7) 

Here, a = 2.9, 2.5, 17.2 for the Plummer, Wq = 2 and 
Wq = 9 King profiles, respectively. As a first-order 
approximation, we can also substitute <j 3 (0)to _ 1 (0) 
from Eq. 2, so that: 

A«/?f rc (0)p 1/2 (0), (8) 


A useful definition of relaxation timescale for a clus¬ 
ter of mass M is, from [14]: 


t r 


4tt(3/2) 1 /2 G 2 (lnA)nm 2 ’ 


(1) 


where [3 = 0.022, 0.015, 0.131 for the three selected 
models, respectively (in CGS units). 

Finally, it can be shown that 


r h (0)- 7 G- 1 / 2 p- 1 / 2 (0) CTl (0), (9) 


where 173 is the three-dimensional velocity dispersion, 
n is the number density of stars, m is the average stel¬ 
lar mass and In A is the Coulomb logarithm. For pro¬ 
cesses related to core collapse, we can take the initial 
values of these quantities at the cluster center (e.g., 
El) ; we obtain a central relaxation timescale t rc : 


f rc (0) 


3 3 / 2 ct?(0) 

47t( 3/2) 1 / 2 G 2 (ln A)p(0)?n(0) ’ 


( 2 ) 


where p( 0 ) is the initial mass density at the core and 
tJi ( 0 ) is the initial one-dimensional velocity dispersion. 
The initial average stellar mass m(0) at the core may 
be larger than the average mass over the whole cluster, 
if there is initial mass segragation. For simplicity, we 
shall assume they are the same. 

Assuming for simplicity that the system is virial- 
ized, we have an additional relation between density 
and velocity dispersion at the cluster center. For the 
Plummer model ( 0 ; 0 ), 

erf (0) = (1/2) 3 / 2 (4tt/3 ) 1/2 G 3/2 Mp 1/2 (0), (3) 


hence 


(7i (0) = 0.518 G 1 / 2 M 1 / 3 p 1 / 6 (0). 

(4) 

Analogous expressions can be derived from the virial 
theorem for the King density profiles ([27|). For ex¬ 
ample, for the King model with dimensionless central 
potential Wo = 2 ll28lb 

(7i (0) = 0.585 G 1 / 2 M 1 /, 3 p 1 / 6 (0), 

(5) 

and 


(7i (0) = 0.286 G 1 /, 2 M 1 / 3 p 1 / 6 (0), 

( 6 ) 


where the constant 7 = 1.56, 1.17, 13.10 for the three 
profiles. 

Using Eqs. (2) and (7) one can now plot the cen¬ 
tral velocity dispersion as a function of central den¬ 
sity, for a given central relaxation timescale: ( 73 ( 0 ) = 
/[p(0);f rc (0)]. Using one of Eqs. (4/5/6), one can 
plot cti(0) = g[p(0);M]. From Eq. 9, one can plot 
(j 1 (0) = h[p( 0); rh(0)]. See also in particular their 
Figure 1.4. 


3. TIMESCALE FOR CLUSTER 
EVOLUTION 

Numerical simulations (0 ) for a variety of Plum¬ 
mer and King profiles have shown that the timescale 
for core collapse (f C c(0)) in clusters with a broad mass 
spectrum is proportional to the central relaxation 
timescale, rather than to the relaxation timescale at 
the half-mass radius: 


t cc (0) « 0.15f rc (0), (10) 

also in agreement with the simulations of [lj. The 
final core mass after the runaway collapse is found to 

be (|T3): 


M cc « 0.002M, (11) 

which is likely to evolve later into a BH with mass 
Mbh ~ 0.001 XI, via direct collapse or SN explosion. 

In the standard treatment, the timescale in Eq. (10) 
is compared with the stellar evolution timescale: core 
collapse and the subsequent formation of an IMBH 
are possible only if f cc (0) 3 Myr, typical lifetime of 
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the most massive stars on the main sequence (fl 2 |). 
If the core collapse is not completed after 3 Myr, it 
will be stopped by strong mass losses from supergiant 
winds and SNe. By imposing this constraint, one can 
identify the regions in the (eri(O), p(0)) plane where 
core collapse is most likely to occur. 

The situation becomes more complicated when the 
role of gas in the young clusters is taken into ac¬ 
count. Both from an analytical approximation and 
from numerical simulations (e.g., |3Cj, in particular 
Figure 1.21) it appears that centrally-concentrated 
clusters can achieve core collapse in < 1 Myr. If so, 
the runaway merger of the most massive stars occurs 
when most stars in the embedded cluster are still sur¬ 
rounded by a spherical cocoon, optically thick to the 
ionizing radiation (class 0/1 protostars). The lifetime 
of a class I phase is « 1-3 x 10 5 yr for low-mass proto¬ 
stars (e.g., 0 ; -12 : : 33 j), and probably even longer 
for O stars, which spend ss 13-15% of their lifetime 
(~ 3-5 x 10 5 yr) shrouded by an optically thick cocoon 
(13; 253 )- Taking into account that even “instanta¬ 
neous” star formation in a young cluster is in fact 
spread out over ~ a few 10 5 yr, a young cluster re¬ 
mains embedded in molecular gas for ss 0.5-1 x 10 6 yr 
(e.g., 0 )- At these early ages, the remaining cold gas 
in the young cluster has a mass at least comparable 
to or larger than the mass in stars. 

When most of the cocoons dissipate, the cluster gas 
is quickly ionized (on a timescale of ~ 10 5 yr) by the 
Lyman continuum photons emitted by the OB stars, 
and reaches a characteristic temperature ss 10 4 K. For 
typical densities ~ 10 3 - 10 6 cm -3 , this corresponds to 
pressures P/fce ~ 10' -10 10 K cm -3 , many orders of 
magnitude larger than the pressure of the ISM. As a 
consequence, the gas in the cluster expands “explo¬ 
sively” with a velocity of order of the sound speed, 
v ~ 10 km s _1 (|22l; 0|)- In some young clusters 
in the Antennae, gas expansion velocities ~ 25-30 km 
s _1 have been inferred (0|; lEoj). 

If the gravitational potential in the cluster is too 
shallow, the expanding gas becomes unbound; if the 
mass loss ^ 50%, the whole cluster will dissipate on 
the same timescale. As a back-of-the-envelope esti¬ 
mate, this occurs when the gas expulsion speed is 
larger than the velocity dispersion at the half-mass 
radius: < 73 (r = rh) 10 km s _1 (0). For typical 
Plummer and King profiles, this corresponds to a cen¬ 
tral velocity dispersion eri(0) ;$ 15-25 km s _1 . A more 
accurate calculation of the exact value of the velocity 
dispersion threshold depends on the details of the gas 
and stellar distribution and is beyond the scope of this 
work. For simplicity, we shall take cri(0) < 20 km s _1 
as a condition for cluster disruption, and (conserva¬ 
tively) t = 5 x 10 5 yr for the lifetime of an embed¬ 
ded cluster. At the same time, clusters have to be 
more massive than ~ 1O 4 M 0 for the explosive disrup¬ 
tion process to occur: smaller clusters do not contain 
enough O stars to ionize all the gas. 


The analytical approximation outlined in Section 2 
did not take into account the role of gas. A simple 
way to account for this effect is to substitute the total 
cluster mass M to t = M + M s = fM for the stellar 
mass M in the virial theorem; / depends on the star 
formation efficiency and cluster age. For the young 
clusters we are dealing with, / ss 1.5-3. As a result, 
(j 1 (0) in Eq. (4/5/6) will be multiplied by a factor 
/% 3 . For example, taking / = 2, we shall plot 

o-i (0) = 0.65 G 1/2 M 1/3 p 1/6 (0) (12) 

for the Plummer model, and 

o-i (0) = 0.36 G 1/2 M 1/3 p 1/6 ( 0), (13) 

for the Wo = 9 King profile, where p(0) = /n(0)m(0) 
is now the total initial central density. With simi¬ 
lar changes we also easily modify the other scaling 
relations in Section 2; for example A = fM/m ~ 
(/?//)M 0 )p 1 / 2 ( 0 ), etc. 

4. IMBHs WITH OR WITHOUT THEIR 
PARENT CLUSTERS 

Putting together the timescale constraints discussed 
in Section 2 and 3, we identify two different regions 
in the parameter space, relevant to the formation of 
IMBHs from runaway core collapse: 

• the first regime (region A in Figs. 1-4) is the one 

discussed by [ 12 : core collapse on a timescale 
t c c (0) 3 Myr, inside a super star cluster (mass 

~ a few 10 5 to a few 1O 6 M 0 ), giving rise to a 
compact remnant 5OOM 0 . The parent cluster 
is massive enough to survive the initial gas ex¬ 
pulsion; it may dissipate later, on much longer 
timescales, owing to mass loss in later stages of 
stellar evolution, or tidal interactions. Stellar 
populations near the brightest ULXs in nearby 
star-forming galaxies have typical ages ~ 10 7 
yr: if the accreting compact objects are IMBHs 
formed through this process, we expect them to 
be still contained inside their parent clusters. A 
likely example is the ULX associated to the su¬ 
per star cluster MGG -11 in M82 ('pH]'): 

• the second regime (region B in Figs. 2-4) is 
for t cc (0) ^ 0.5 Myr, in clusters that satisfy 
the disruption condition discussed in Section 3: 
M > 1O 4 M 0 and <ti( 0) < 20 km s _1 . A cluster 
in this region of the parameter space is unlikely 
to survive beyond its embedded phase: its stars 
will keep expanding freely with velocities ~ ( 73 . 
After 10 7 yr, this will result in an OB associ¬ 
ation with a diameter of ~ 200 pc, difficult to 
distinguish from other surrounding star-forming 
regions. Not all the clusters in this regime are 
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suitable for IMBH formation: assuming that 
Eq. (11) holds, only the subset of exploding clus¬ 
ters with M 5 x 10 4 Mq are massive enough 
to produce BHs with M 50M 0t required to 
explain the observed ULX luminosities. 

The densest young star clusters in our galaxy (e.g., 
RU136, NGC 3606, the Arches) have total central den¬ 
sities ~ 1O 6 M 0 pc -3 . Therefore, we took p max (0) = 
1O 7 M 0 pc -3 as a reasonable upper limit to the initial 
central density in our parameter space. 

Core collapse is more likely to occur on short 
timescales for centrally concentrated profiles: for ex¬ 
ample, for the Wo = 9 King profile (Fig. 2) rather than 
for the Plummer profile (Fig. 1), in agreement with 
the simulations of |3Cj. For a given cluster model, the 
parameter space available for runaway core collapse 
is enhanced when the average stellar mass is higher, 
that is when the initial mass function is top-heavy. 
Here we compare the initial conditions for m = 0.5M 0 
(Fig. 3) and m = 3.0M 0 (Fig. 4), for the same Wo = 9 
King profile. For example, values of m ~ 3.0M 0 have 
been inferred for the MGG-11 cluster in M82, which 
is likely to contain an IMBH ([b|; |uj). 

Taking a W = 9 King profile with a typical central 
density ~ 1O 6 M 0 pc -3 (Figs. 2-4), we infer from our 
simple analytical approach that clusters with initial 
stellar masses ~ 1O 5 M 0 may produce a collapsed core 
during their embedded phase, but may not survive the 
ionization/expulsion of their gas. The collapsed cores 
may later evolve into apparently isolated IMBHs with 
masses up to « 2OOM 0 . On the other hand, super 
star clusters with initial stellar masses ~ 1O 6 M 0 may 
survive the embedded phase and produce IMBHs with 
masses up to ~ 1O 3 M 0 . Considering that the initial 
cluster mass function diV ~ M~ 2 dM, we also expect 
that IMBHs from clusters in region B should be an 
order of magnitude more numerous than those formed 
from clusters in region A. 

5. CONCLUSIONS 

We have used a simple analytical approximation to 
suggest that runaway core collapse in a young star 
cluster may occur in two distinct regimes. The first 
regime (extensively investigated with numerical sim¬ 
ulations by 0 and 0 ) is a collapse on a timescale 
;$ 3 Myr (main-sequence lifetime of the O stars) in a 
cluster of total initial mass ~ 10 6 M 0 . Observed after 
~ 10 7 yr, the outcome of the collapse is likely to be 
an accreting IMBH (Mbh ~ 1O 3 M 0 ) inside a bright 
cluster. The second regime is a collapse on a shorter 
timescale (^ 0.5 Myr), in a smaller cluster (initial 
mass ~ 10 5 M 0 ), during its initial embedded phase, 
when the most massive protostars are still shrouded 
by optically-thick cocoons. If the velocity dispersion 
in a cluster is smaller than the thermal velocity of the 


ionized gas, such a cluster does not survive beyond its 
embedded phase, because of the explosive loss of gas. 
Observing the system after ~ 10 7 yr, we may find an 
accreting IMBH (Mbh ~ 1O 2 M 0 ) in a star-forming 
region or OB association, but apparently not inside 
any clusters. 

Observationally, ULXs have been explained as ac¬ 
creting IMBHs. Runaway core collapse is a possible 
mechanism of IMBH formation in young stellar envi¬ 
ronments. However, most ULXs are not coincident 
with a cluster, though they are often associated with 
OB stars. We have suggested that the existence of 
two possible regimes of core collapse may explain this 
puzzle. The ULX associated with the M82 cluster 
MGG-11 may be an example of core collapse in the 
first regime, when the parent cluster survives. Most of 
the other ULXs in nearby star-forming galaxies, with 
an OB companion and X-ray luminosities up to 10 4 ° 
erg s , may have been formed in the other regime, 
when the parent cluster evaporates explosively. 

The parameter space allowing IMBH formation 
from a runaway core collapse is more extended for 
a top-heavy stellar mass function. It has often been 
claimed that the stellar initial mass function is top- 
heavy (higher fraction of high-mass stars) in some 
starburst galaxies, although the issue is still contro¬ 
versial (|42} and references therein). If this is the case, 
we speculate that it may explain why ULXs seem to 
be more often found in starburst galaxies or in envi¬ 
ronments where star formation is triggered by galactic 
mergers or tidal interactions. 

Finally, if many ULXs were indeed formed in the 
core of a long-dissolved protocluster, we may want to 
search for clues of that disruption event. For example, 
after ~ 10 7 yr, the expanding shell or cloud of gas 
with a mass ~ 1O 5 M 0 will have a characteristic radius 
~ 100 pc and density ~ 1 cm -3 . This gas may now 
be X-ray photoionized by the accreting IMBH. We 
speculate that at least some of the ionized nebulae 
discovered around nearby ULXs (@|) may be related 
to the disruption of the same protocluster where the 
BH was formed. 
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Figure 1: Parameter space for the initial central density p( 0) and one-dimensional velocity dispersion <ti( 0), in a cluster 
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t cc (0) = 0.15f r c(0). A red box (marked with “A”) identifies the region of the parameter space where the runaway core 
collapse occurs in < 3 Myr and the cluster is massive enough to produce an IMBH (Mbh ~ 10 -3 M <1 50Mq). Finally, 
we have divided the parameter space into two subsets with <ti(0) > 20 and < 20 km s _1 . Clusters above the threshold 
are more likely to survive their embedded phase; clusters below the threshold (but more massive than ~ 10 4 Mq) are 
more likely to evaporate, turning into OB associations. Finally, we have chosen an initial density p(0) = 10' Mq pc~ 3 
as a plausible upper limit, based on the values estimated for the densest clusters in our and nearby galaxies. 
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Figure 2: As in Figure 1, for a King profile of concentration index Wo = 9. Here there are two distinct regions of the 
initial parameter space where runaway core collapse can occur: region A is for massive clusters which are expected to 
survive bound, and region B for less massive systems which are expected to evaporate. We speculate that IMBHs 
formed from region-B clusters may explain a population of ULXs with X-ray luminosities L x « 10 40 erg s~ 4 found in 
nearby star-forming galaxies but not associated to any present-day cluster. 
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Figure 3: As in Figure 2, for an average initial stellar mass m = 0.5Mq. The initial parameter space available for 
runaway core collapse and IMBH formation is much reduced. 
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Figure 4: As in Figure 2, for an average initial stellar mass m = 3.0 Mq. The initial parameter space available for 
runaway core collapse and IMBH formation is enhanced for a top-heavy stellar mass function. 












